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INTRODUCTION 


A missile  can  kill  a target  in  several  ways,  for  example,  by  an 
internal  or  external  blast  from  the  warhead  and  by  fragments  from  the 
warhead.1  A mathematical  model  for  determining  the  single-shot  kill 
probability  requires  that  the  path  of  the  fragments,  the  shape  of  the 
blast  envelope,  the  shape  of  the  target,  etc.  be  expressed  mathemati- 
cally. J.  von  Neumann  suggested  that  a set  of  ellipsoids  could  be  used 
to  define  the  blast  envelope  and  the  shape  of  the  target.2  For  swept- 
back  configurations  of  the  target,  a set  of  cardiaci  can  be  used  to 
approximate  the  target  mathematically.3  The  initiation  of  certain 
fuzes  on  the  missile  require  knowledge  of  the  distance  between  the 
target  and  the  missile's  fuze.  Mathematically,  this  is  the  distance 
between  the  fuze  and  the  nearest  ellipsoid  describing  the  target. 


METHOD  OF  THE  CALCULUS  OF  VARIATIONS 


We  wish  to  find  the  shortest  distance,  d,  between  a fixed  point, 
P(x<>,yo,Zo).  and  the  surface,  G(x,y,z)  = 0.  Q(x,  ,y,  ,Zi)  is  defined  to 
be  a point  on  this  surface  such  that  the  distance  from  P to  Q is  a 
minimum.  See  Figure  1. 


Figure  1.  Problem  Configuration 
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o (5  the.  Single-Shot  Kill  Probabilities  o<5  a General  Missile  Versus 
a General  Aircraft,"  Ballistic  Res ear.ch  Laboratories  Memorandum 
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3.  V.  Kucher,  " Mathematical  Approximation  oh  Sweptback  Conhlguratlons 
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We  shall  consider  the  surface  to  be  an  ellipsoid  such  that 

2 2 2 

t>(x,y , z)  = — + £-  + — -1  = 0,  a,b,c  > 0.  (1) 

a2  b2  c2 

Any  ellipsoid  that  is  not  expressed  in  standard  form  can,  by 
translation  and  rotation  transformations,  be  written  in  the  form  of 
Equation  1.  Naturally,  these  transformations  would  also  be  applied 
to  the  fixed  point  in  space. 

The  length  of  any  curve  extending  from  point  P to  point  Q is 
given  as 


S = 


rx  i 


[i  + Cy ' ) 2 + (z')2]  2 dx. 


xo 


A procedure  for  minimizing  an  integral  of  the  form4 


fxl 

1 = F(x,y,z,y' ,z')dx 


is  to  apply  the  Euler-Lagrange  equations: 


IE  . iLfiE  ) = o 

3y  d?3y'J  u 


and 


(2) 


(3) 


(4) 


3F 

3z 


d,  3F 
ebe  3z 


,) 


= 0. 


(5) 


4. 


R.  CouAant  and  V.  HUbvrt,  MeXhodi  ok  Mathematical  PhuUu. 
Mol.  I,  New  Voik,  InteMciience.  PubtUkeM,  1953. 


From  Equation  (2)  we  see  that  for  our  case 


The  Euler-Lagrange  equations  take  the  form 


The  solution  to  the  simultaneous  equations  [ (7)  and  (8)]  is 


where  Cj,  C2,  C3,  and  C4  are  constants  of  integration.  As  expected 
we  have  found  that  the  shortest  curve  between  points  P and  Q is  a 
straight  line. 


Applying  the  conditions  at  point  P that  x = xq,  y * yo»  and 
:n,  we  have 


The  next  step  in  our  procedure  to  minimize  the  distance  between 
points  P and  Q is  to  apply  the  transversality  conditions: 


The  slope  of  the  line  given  by  Equations  (9)  and  (10)  is  provided 
by  Equations  (19)  and  (20) . Thus 


a2y  a2z 

Cx  = - and  C3  = (21)  - (22) 

b2xj  c2Xj 


Note  that  the  slope  of  the  surface  of  the  ellipsoid  at  point  Q can  be 
determined  to  be  the  negative  reciprocal  of  the  slope  given  by  Equations 
(19)  and  (20);  consequently,  the  line  is  perpendicular  to  the  surface  of 
the  ellipsoid.  Thus  transversality  coincides  with  orthogonality. 


Substituting  Equations  (21)  and  (22)  into  Equations  (11)  and  (12), 
we  find  that 


a 2V\ 


c2  = yo  — - — *o» 

b2x. 


(23) 


a*-z. 


C4  = zq  - 


c2x. 


x0. 


(24) 


Substituting  the  four  constants  of  integration  into  Equations  (9) 
and  (10),  we  obtain  the  equation  of  the  shortest  curve  between  points 
P and  Q: 

a2y  j 

y = (x  - x0)  + y0 , (25) 

b2Xj 


a l j 

z = (x  - x0)  + Z0.  (26) 

r^v . 
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Evaluating  Equations  (25)  and  (26)  at  Q,  we  find  that 


b2yoxl 


yi 


zi 


b2xj  - a2(xj-x0) 


C2Z0X! 


c2xx  - a2(Xi-x0) 


(27) 


(28) 


If  we  evaluate  Equation  (1)  at  Q and  substitute  Equations  (27)  and 
(28)  into  the  resulting  equation,  we  eliminate  y!  and  zj  and  obtain  a 
sixth-degree  polynomial  in  xj: 


K0xi  + K^xi  + K2X1  + K3X1  + K4X1  + K5X1  + K6  = 0 (29) 


with  coefficients  given  by 


K0  = r2s2. 


Kj  = 2rskt, 

K2  = k2u  + a2(gs2  + hr2  - r2s2) , 
K3  = 2k( k2t  + a2(gs  + hr  - rst)] , 
K4  = k2[  k2  + a2(g  + h - u)] , 

K5  = -2a2k3t, 

K6  = -a2k1*. 
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The  reader  may  use  his  favorite  method  or  computer  routine  to  find 
the  roots  of  Equation  (29).  Since  Kg  is  negative  and  Kq  is  positive, 
we  are  assured  that  at  least  one  positive  real  root  and  one  negative 
real  root  exist.  The  fact  that  IxJ  < a is  useful  in  the  determination 
of  the  roots  which  have  meaning  in  this  problem. 

If  point  P(xo,yo»zo)  does  not  lie  in  any  of  the  coordinate  planes, 
values  of  yj  andzi,  corresponding  to  each  of  the  real  roots  (Ixj]  < a) 
of  Equation  (29),  are  obtained  from  Equations  (27)  and  (28).  Point 
Q (X1 , y l , z i ) is  found  by  determining  the  minimum  value  of  d for  each 
triplet  (xi,yi,zi)  obtained,  where 

d = [ (Xj  - x0)2  + (yj  - y0)2  + (zx  - zq)2]^.  (30) 


As  an  example  of  the  above  procedure,  let  (a,b,c)  = (4,3,2)  and 
(x0>yo>zo)  = (1»1»1).  The  roots  of  Equation  (29)  for  this  example 
were  determined  by  an  available  computer  routine  based  on  the  procedure 
developed  by  Bareiss  and  Hamelick.5  However,  other  techniques,  such  as 
the  Newton-Raphson  method,  should  be  adequate  for  this  purpose.  A 
summary  of  the  calculations  is  listed: 


5.  E.  8 aAeiAA  and  R.  Hamelink,  "RSSR  Routine.,  A Root-SquaAing  and 
SubAeAultant  PAOceduAz  fioA  Rinding  Zeno  A oi  R exit  Polynomiath , " 
AAgonne  National  LaboAatoAy  RepoAt  No.  ANL-69&7,  OctobeA  1 965. 
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Root 

1 

2 


y i 


*1 


1.118,521,8  1.232,101,4  1.735,663,8  0.780,461,2  (minimum) 

-3.821,653,6  -0.804,529,0  -0.247,116,9  5.297,166,1 


The  procedure,  which  was  developed  for  determining  Q and  d,  may 
require  some  modification  if  point  P lies  in  any  of  the  coordinate 
planes.  For  example,  let  P be  a point  such  that  Yq  = 0 and  Xo,zo  / 0. 
Equation  (27)  would  force  y^  to  be  zero.  With  z\  determined  by 
Equation  (28),  it  is  possible  that  the  point  (x!,0,zi)  will  not  lie  on 
the  surface  of  the  ellipsoid. 

Rather  than  modify  the  above  procedure  for  cases  where  P lies  in 
any  of  the  coordinate  planes,  the  problem  will  be  analyzed  by  the 
method  of  Lagrange  multipliers. 


METHOD  OF  LAGRANGE  MULTIPLIERS 

Let  the  square  of  the  distance  between  P and  any  point  (x,y,z)  on 
G,  Equation  (1),  be  given  as 


D = (x  - x0)z  + (y  - 
Define  the  function,  f,  as 

f = D 

where  X is  a Lagrange  multiplier. 

The  solutions  of  the  equations 

3f 
3x  ' 

3f  _ 

3y  = 

af 

3z  = 

are,  respectively, 


y0) 2 * (z  - z0)z.  (31) 

XG,  (32) 

0, 

0,  (33) 

0, 
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x(l  - — ) = x0, 
a2 

y Ci  - — ) = yo» 

b2 


Z(1  - 


Zq  . 


(34) 


The  set  of  Equations  (34)  and  Equations  (1)  and  (30)  are  used 
for  determining  Q and  d when  P is  in  any  of  the  coordinate  planes,  as 
illustrated  in  the  example  that  follows. 


Equations  (34)  take  the  form 


* 0. 

Let 

a,  b 

x(l  - 

= *0 

a2 

y ( i - 

b2 

= o. 

z(l  - 

= 0. 

If  A ^ b2,  c2.  Equations  (35)  yield 


(35) 


(36) 


Equations  (1)  and  (36)  yield 

A = a2  + ax0 

and  thus  x = ±a.  Consequently,  a possible  choice  for  Q is  (±a,0,0). 
If  A = b2  t c2.  Equations  (35)  yield 

z = 0, 

a2x0 

x = . 

a2  - b2 
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If)  i win 


- . ....  . ...... 


In  summary,  if  P is  (xo,0,0),  x0  # 0,  and  a,  b,  and  c are  distinct, 
the  candidates  for  Q are  the  following: 

1.  (±a,0,0). 


Point  Q is  selected  by  the  substitution  of  the  above  coordinates 
into  Equation  (30)  and  by  examining  the  resulting  values  of  d for  a 
minimum . 

Working  with  the  basic  Equations  (34) , we  may  determine  possible 
candidates  for  Q when  P is  in  any  of  the  coordinate  planes.  A summary 
of  possible  coordinates  for  Q for  various  positions  of  P in  the  coordi- 
nate planes  are  given  in  Table  I in  Appendix  I. 

If  several  sets  of  coordinates  for  Q are  listed  for  a particular 
case  in  Table  I,  Equation  (30)  is  used  to  determine  which  set  yields 
the  minimum  value  for  d;  however,  several  cases  exist  where  there 
could  be  an  infinite  number  of  possible  positions  for  Q.  Case  2, 
for  example,  shows  that  Q could  lie  on  a circle  located  in  the  x,y- 
plane.  By  selecting  any  point  on  this  circle,  we  can  calculate  d. 

Several  cases,  such  as  Case  21,  require  that  roots  of  quartics 
be  determined.  Again,  the  reader  has  the  option  to  use  his  favorite 
scheme  to  obtain  these  roots. 


OTHER  METHODS 

The  dynamic  programming  technique  and  the  method  of  Lagrange 
multipliers  were  used  to  find  approximate  solutions  to  the  problem  of 
the  distance  between  a point  and  an  ellipsoid6.  The  dynamic  programming 
approach  considered  the  axes  of  the  ellipsoid  to  be  parallel  to  the 
coordinate  axes;  the  Lagrange  multiplier  approach  considered  the  case 
where  these  axes  did  not  necessarily  have  to  be  parallel. 

It  is  not  clear  from  the  dynamic  programming  solution  as  to  how 
situations  such  as,  for  example.  Cases  2 to  5 in  Appendix  I,  where  an 

6.  A.  CeZmint  and  W.  Sacco,  "NumeAlcat  Computation t o&  the  Pittance 
h*jom  a Point  to  an  EUlptold,  " SaJUUttic  ReteaAch  LabotatosUet 
Report  No.  132S,  1966.  (AD  #641011) 
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infinite  number  of  points,  Q,’  solve  the  problem,  are  resolved.  The 
computer  algorithm,  developed  from  the  Lagrange  multiplier  method, 
considered  P to  be  outside  the  ellipsoid;  otherwise,  the  algorithm 
became  too  complicated. 


This  report  presents  exact  solutions  when  possible;  otherwise, 
the  solutions  are  obtained  by  finding  the  roots  of  6-th  or  4-th  degree 
polynomials  in  a known  domain  for  the  roots.  For  a given  accuracy,  a 
computer  solution  of  the  approach  presented  in  this  report  should  be 
faster  than  those  in  Reference  6.  Furthermore,  since  P can  be  located 
inside  or  outside  of  the  ellipsoid,  all  requirements  for  the  application 
of  this  study  to  systems  analysis  have  been  satisfied. 


SUMMARY 


We  wish  to  determine  the  shortest  distance  between  a fixed  point 
P(xo»yo»zo)  and  an  ellipsoid.  Q(xi,yi,zi)  is  defined  to  be  on  this 
surface  such  that  the  distance  between  P and  Q is  the  minimum. 

Translation  and  rotation  transformations  should  be  used  on  the  system 
consisting  of  a fixed  point  and  an  ellipsoid  to  describe  the  ellipsoid  in 
standard  forfo,  Equation  (1),  The  real  roots  of  Equation  (29),  such  that 
|xj|  < a,  with  their  corresponding  components,  given  by  Equations  (27) 
and  (28) , provide  possible  values  for  Q if  P does  not  lie  in  any  of  the 
coordinate  planes.  If  P does  lie  in  one  or  more  coordinate  planes. 
Appendix  I provides  possible  coordinates  for  Q.  In  either  case,  the 
final  choice  of  coordinates  for  Q are  determined  by  finding  the  minimum 
value  for  d by  Equation  (30) . 
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